Finite Volume Chiral Partition Functions and the Replica Method 



Jesper Levinseqj 
The Niels Bohr Institute 
Blegdamsvej 17, DK-2100 Copenhagen 
Denmark 
(Dated: February 1, 2008) 

In the framework of chiral perturbation theory we demonstrate the equivalence of the super- 
symmetric and the replica methods in the symmetry breaking classes of Dyson indices j3 = 1 and 
(3 — 4. Schwinger-Dyson equations are used to derive a universal differential equation for the finite 
volume partition function in sectors of fixed topological charge, v. All dependence on the symmetry 
breaking class enters through the Dyson index /3. We utilize this differential equation to obtain 
Virasoro constraints in the small mass expansion for all (3 and in the large mass expansion for /3 = 2 
with arbitrary v. Using quenched chiral perturbation theory we calculate the first finite volume 
correction to the chiral condensate demonstrating how, for all /3there exists a region in which the 
two expansion schemes of quenched finite volume chiral perturbation theory overlap. 
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I. INTRODUCTION 

In the low-energy and chiral limit of QCD, chiral sym- 
metry is spontaneously broken. The physics is dominated 
by the pseudo-Goldstone bosons, and the corresponding 
effective field theory is known as chiral perturbation the- 
ory. This theory is describable in terms of a chiral La- 
grangian. However, the physics in this low-energy regime 
depends upon the exact way in which the chiral symme- 
try is spontaneously broken. There are believed to be 
three ways in which spontaneous chiral symmetry break- 
ing can happen []}. These depend upon the representa- 
tion of the fermions in the following way: 

• The representation of the fermions is pseudo- 
real. The expected symmetry breaking pattern is 
SU(2N f ) - Sp(2N f ). 

• The representation is complex. In this case 
we expect the symmetry breaking pattern to be 
SU L (N f ) x SU R (N f ) SU(N f ). An example 
is fermions in the fundamental representation with 
the number of colours, N c > 3. Thus ordinary QCD 
falls into this class. 

• The fermions are in a real representation. The 
expected symmetry breaking pattern is then 
SU(N f ) - SO(N f ). 

These symmetry breaking patterns are nowadays labeled 
by the Dyson- indices, f3 = 1, (3 = 2, and /3 = 4, respec- 
tively, due to a connection to Random Matrix theory 
0. For a symmetry breaking pattern G — > H, the fields 
live on the coset G/H. The usual perturbation scheme 
for chiral perturbation theory is for large volumes and is 
an expansion in terms of the momenta of the Goldstone 
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modes. This perturbative scheme will be applied in most 
of this paper. 

A systematic approach to calculations in (partially) 
quenched chiral perturbation theory is known as the su- 
persymmetric method 0, Gl E|. I n t ne supersymmet- 
ric formulation k "valence" quark species are introduced 
with k ghost quarks of opposite (bosonic) statistics in 
addition to the Nf "sea" quarks. In this way the effec- 
tive partition function is extended, becoming a genera- 
tor of n-point functions with the additional quarks act- 
ing as source terms. The terminology "supersymmetric 
method" reflects the fact that the chiral flavor symmetry 
group is extended to a super Lie group, while space-time 
supersymmctry is not intrinsic to the method. 

The replica method, which shall be applied in this pa- 
per, has turned out to be an alternative to the supersym- 
metric method. In this method one adds N v fermionic 
valence quarks taking N v to zero at the end of the cal- 
culations. Obviously, if we let N v — ► we recover the 
original partition function. However, as in the super- 
symmetric method this extension of the partition func- 
tion makes it a generating functional for n-point func- 
tions with the sources being the valence quarks. It can 
be convenient to add k sets of N v valence quarks with 
masses m Vl , . . . , m Vk , but for the applications in this pa- 
per we need only include one set of N v valence quarks of 
mass m v . We assume the symmetry breaking patterns 
in this extended theory to be the usual symmetry break- 
ing patterns. Thus for intermediate calculations in this 
paper, in the above-mentioned symmetry breaking pat- 
terns we will simply replace Nf by Nf + N v . Since the 
limit N v —* formally requires an analytic continuation 
of the Lie groups to non-integer dimensions the validity 
of the method is not obvious, but it can be shown that 
the required analytic continuation can be performed in 
series expansions |H S @ • Recently there have been some 
developments into non-perturbative results in this frame- 
work [3, but these are out of the scope of this paper. In 
both this framework and the supersymmetric one, the 
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fully quenched limit corresponds to vanishing Nf while 
the theory is partially quenched for Nf non-vanishing. 

The equivalence of these methods have been demon- 
strated for (3 = 2 and thus for ordinary QCD |(| by ex- 
plicitly showing the equivalence of the Feynman rules at 
the one-loop level. In the first part of this paper we will 
extend this analysis to the other classes of chiral symme- 
try breaking, showing the equivalence only explicitly at 
the one-loop level but noting that the analysis is easily 
carried over to the higher loop level. To illustrate the 
equivalence we will also calculate the partially quenched 
chiral condensate at one-loop order, illustrating how the 
statistics signs of the supersymmetric formulation are re- 
produced in the replica formulation by the counting of 
possible quark-loops. 

In the second part of this paper we present one of the 
main results of this paper, a universal second order differ- 
ential equation to determine the finite volume chiral par- 
tition function in all classes of chiral symmetry breaking 
and for arbitrary topological charge, v. This differential 
equation is written in terms of the fermion masses and 
all dependence on the pattern of symmetry breaking is 
through the Dyson index, (3. The form of the differential 
equation in the specific case of (3 = 2 with zero topologi- 
cal charge is already known |l(lj | , but the generalization of 
this paper has not been found previously. The method 
is to develop Schwinger-Dyson equations for the effec- 
tive low energy, or finite volume, partition function and 
next to realize how these can be written in terms of the 
fermion masses. For (3 = 2 with vanishing topological 
charge the effective partition function has been recog- 
nized as a group integral of Kontsevich type with 
potential V{X) = 1/X and this model, also known as 
the one-link integral, has been intensively studied in the 
literature, see, e.g., refs. [1 111 HI llH . 

Having found a governing differential equation we pro- 
ceed to solve this equation. It should be realized that 
there are two relevant expansions, the small-mass and 
the large-mass expansions. Using the invariance prop- 
erties of the effective finite volume partition function it 
is possible to find Virasoro constraints. For the small 
mass phase this was first done in ref. 1 1 11 ] in the sim- 
plest version, namely (3 — 2 with vanishing topological 
charge, by recognizing the partition function as being a 
unitary integral of Kontsevich type. More recently, Vi- 
rasoro constraints in the small mass phase for general 8 
and topological charge have been calculated in ref. Q 
by using a method much similar to the one presented 
here. We exactly reproduce the Virasoro constraints of 
ref. 0- The large mass expansion is more involved as 
it is not possible to do a simple perturbative expansion, 
instead the large mass expansion is an asymptotic, or 
saddle-point, expansion 10, 11]. The saddle-point ap- 
proximation corresponds to the classical limit and is an 
expansion in powers of l/Nf. In the large mass limit 
Virasoro constraints have been previously determined in 
the case of (3 = 2 0] . We find Virasoro constraints in 
the general case of j3 = 2, where the inclusion of a non- 



vanishing topological charge is thus a new result. At next 
to leading order in l/Nf the saddle-point approximation 
turns out to be a simple correction to the leading order 
term in addition to an (apparently) infinite expansion in 
the inverse fermion masses. However, this last expansion 
is proportional to 1 — 2/ j3 and thus it vanishes if (3 = 2. 
Precisely this property allows us to calculate large-mass 
Virasoro constraints in the (3 = 2 case with arbitrary 
topological charge. Although we do not determine Vira- 
soro constraints in the other classes of symmetry breaking 
we find that it is nevertheless possible to extract useful 
information from the saddle-point approximation; in the 
case of equal fermion masses the differential equation for 
the partition function is still tractable. We solve the gov- 
erning differential equation for equal fermion masses and 
as an application we find the lowest order corrections to 
the fully quenched chiral condensate in sectors of topo- 
logical charge, a result which will be very useful in the 
last part of this paper. 

As Gasser and Leutwyler 0] were the first to point 
out, there are two finite volume regimes to consider. The 
first is the large, but finite, volume which can be con- 
sidered a small perturbation to the infinite volume the- 
ory. This is the volume considered in the first parts 
of this paper and to which the usual chiral perturba- 
tion theory in terms of a momentum expansion of the 
chiral Lagrangian applies. The last part of the paper 
concerns quenched chiral perturbation theory in volumes 
much smaller than the correlation lengths of the Gold- 
stone modes. Although this may seem ill-defined at first 
it is still possible by means of numerical simulations to 
extract physical, infinite volume quantities from this the- 
ory. In this volume-regime the usual momentum expan- 
sion breaks down due to the presence of zero-momentum 
modes 16]. The different expansion scheme required in 
this regime, known as the e-expansion |15| . treats this 
problem by means of a collective field technique, col- 
lecting the zero-momentum modes in one field and the 
non-zero modes in another. The non-zero modes are still 
perturbative. Difficulties in this scheme arise from the 
fact that the integration over the zero-momentum modes 
has to be done exactly. 

In ref. ^ was shown that in the case of (3 = 2 it is 
possible to determine a region in which the two pertur- 
bative expansion schemes overlap also in the quenched 
limit. We extend this analysis to the other classes of 
chiral symmetry breaking, seeing how also in these cases 
such a region exists and demonstrating how the quenched 
chiral condensates of the two regimes exactly match in 
the above-mentioned region. 

This paper is organized as follows. In Section [H] we 
consider (partially) quenched chiral perturbation theory. 
The equivalence of the supersymmetric and the replica 
method in the alternative ((3 = 1,4) classes of spon- 
taneous chiral symmetry breaking is explicitly demon- 
strated at the one-loop level. Since equivalence was 
demonstrated for ordinary QCD ((3 = 2) in ref. [g, these 
two approaches to calculations in (partially) quenched 



chiral perturbation theory are equivalent in all classes of 
spontaneous chiral symmetry breaking. 

Section II I II concerns finite volume chiral perturbation 
theory. Schwinger-Dyson equations are derived for the 
finite volume chiral partition functions and utilized to 
obtain Virasoro constraints in the limit of small quark 
masses and also (for f3 = 2) in the limit of large quark 
masses. For general (3 we find the next to leading order 
term in a mass expansion of the partition function. Al- 
though the calculations of this section are not considered 
to be in a quenched setting, the relevance of these results 
in such a setting is demonstrated by calculating the first 
correction in a mass expansion to the quenched chiral 
condensate. 

We return to considering quenched chiral perturba- 
tion theory in section IIVI Here it is demonstrated 
how a volume-regime exists in which the two perturba- 
tion schemes of finite volume chiral perturbation theory 
match in the quenched approximation. 



II. THE EQUIVALENCE OF THE 
SUPERSYMMETRIC AND THE REPLICA 
METHODS 

It has been shown that an alternative approach to- 
wards calculations in partially quenched chiral pertur- 
bation theory is offered by the replica method. The 
equivalence between the supersymmetric method and the 
replica method has been demonstrated for /3 = 2 in 
ref. In chiral perturbation theory the natural way 
of demonstrating the equivalence between the methods 
is demonstrating how the generating functional in the 
replica formulation results in Feynman rules equivalent to 
those obtained from the generating functional in the su- 
persymmetric formulation. In this section we will extend 
this equivalence to the other classes of chiral symmetry 
breaking. As was the case in ref. [|| we will only demon- 
strate this equivalence at the one-loop level and thus the 
0(p 2 ) expansion of the Lagrangian density is sufficient 
with the 0(p ) terms acting as counter terms. To demon- 
strate the equivalence, we compare our results with those 
obtained in ref. [lsf by means of the supersymmetric 
method. The extension to the 0(p 4 ) Lagrangian is not 
very difficult, as it only affects the vertices, and not the 
propagator structure. 

In the framework of chiral perturbation theory, the ef- 
fective partition function is given by 



5 -- I dU e -I d4 ^« 

'ueG/H 



(1) 



with the usual effective chiral Lagrangian density to sec- 
ond order in momenta 

ArfF = ^Md^d,U)-^TT(M^(U + U^)) 



The Goldstone fields U are parameterized as 

U = exp(i\/2$/ J F». (3) 

The last two terms in (|2J) contain the flavor singlet field 
$o = Tr$. In the supersymmetric formulation this term 
is an invariant of the reduced symmetry group and it has 
to be included in the replica Lagrangian as well, since 
otherwise the quenched replica limit would not exist |27| . 
For the time being we will ignore the term proportional 
to a and reinstate it at the end of this section by substi- 
tuting TOq with TOq + ap 2 . 

The pion decay constant Fp is chosen so that the fields 
satisfy the usual Gell-Mann-Oakes-Renner relation for 
the masses of the mesons 



Ml = (m l + mj )^. 



(4) 



With our choice of normalization of the fields it turns out 
that 
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Fi = F/y/2. 



(5) 
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Notice that the fact that our choices of normalization of 
the fields and of Fp differ from those of ref. does 
not stem from any peculiarities of the replica method. 
Rather, we find these choices more natural, since they 
ensure that all complex fields are normalized in the same 
way. 

The infinite volume chiral condensate is denoted by 
So and in the replica formulation the partially quenched 
chiral condensate of the additional fermionic copies is, to 
one-loop order, given by 
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= 1 
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The expected symmetry breaking of fermions in 
a pseudo-real representation is SU(2Nj + 2N V ) — > 
Sp(2N f + 2N V ) with Dyson index (3 = 1. 

Following ref. Q] we choose an explicit representation 
of the field $ as 
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1 / v 

71 w ^ 
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The normalization is taken for convenience. In this rep- 
resentation the field (f> is Hcrmitian while the field ip is 
complex and anti-symmetric. Both fields are written in 
terms of (Nf + N v ) x (Nf + N v ) matrices. 
The mass matrix can be chosen diagonal with 



M {1) =\ diag( 



mi, 
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where this normalization ensures that the Gell-Mann- 
Oakes-Renner relation 10} is satisfied. 

To calculate the chiral condensate to one-loop order 
we evaluate the Lagrangian to second order in the fields. 
Only the term containing m 2 , mixes states and this only 
applies to the diagonal fields. This makes it easy to read 
off the propagators for the "off-diagonal" mesons. These 
are seen to be 
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The diagonal mesons also inherit this part, so that their 
propagator becomes 

= ^(p 2 +M 2 ) + 2m 2 /N c . (11) 

This expression can be inverted by noting that for general 
matrices A,B, with Aij = a constant, and Bij = Sijbu 
(no sum over i) diagonal, 

(A + B)- 1 = B- 1 -B~ 1 AB- \ ^ — . (12) 



Thus we find that 

Gi 

where 
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since the masses of the mesons in the replica set are iden- 
tically m v . These expressions for the propagators can be 
compared with the corresponding ex pres sions from the 
supersymmetric method obtained in 18]. For N v — > 
the expressions are equivalent, the difference arising from 
our alternative normalization conventions, only. 

We now proceed to calculate the chiral condensate to 
one-loop order. It is given by J7J) and can thus be cal- 
culated directly from the second order expansion of the 
Lagrangian. We find that 
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where the pion propagator is given by 
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and where we in the last step for simplicity have put 
the sea quark masses equal. This result is identical to 
the result obtained from the supersymmetric formulation 





B. f3 = A 

For fermions in a real representation the expected sym- 
metry breaking class is SU{N f + N v ) -> SO(N f + N v ) 
corresponding to the Dyson index (3 — 4. In this case an 
explicit representation of the Goldstone bosons is ^^] 

/ An -7= A 12 ■ ■ ■ -^Ai^ Nf+Nv \ 



$ = 



\/2 



A, 



.(17) 



Wf+N v ,N f +N v J 



$ is a real, symmetric matrix 0. Thus the normaliza- 
tion of the off-diagonal mesons differs from that of the 
diagonal mesons to ensure correctly normalized kinetic 
terms in the Lagrangian. The mass matrix is diagonal 
with 



= diag(mi, 



,m Nf ,m v , 



(18) 



Following the steps above we evaluate the Lagrangian 
(0) to second order in the fields to obtain the one-loop 
correction to the chiral condensate. As for (3 = 1 the 
propagators of the off-diagonal mesons are easy to read 
off. Again we find 
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The propagators of the diagonal mesons are 

Grl = Satf + M%)+ ml/ N c 
which, inverted by the use of eq. (|12fl . yields 



Gij — 



m 2 /N c 



P 



M 2 



(p 2 + M 2 )(p 2 + M 2 ,)^= 4 (p 2 ) 



with 
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P 



M 2 
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Ml k 
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In the limit N v — > this result also agrees with the result 
obtained in ref. [Tsj from the supersymmetric formula- 
tion. 
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The calculation of the partially quenched chiral con- 
densate proceeds as before. From Q we find 



T,{m v ,{m}) 1 1 at a 
E = l ~F* \2 NfK 



-IJ2 G w(p 2 )) , (23) 
p / 



with the pion propagators given as in eq. (|16() and the sea 
quark masses equal. Again, this result exactly matches 
the result from the supersymmetric method |l8j |. 



C. The chiral condensate in all classes of chiral 
symmetry breaking 

Including the chiral condensate in theories with 
fermions in a complex representation ((3 = 2) [T^. |20|. 
allows us to write, for general (3, 



T,^(m v , {m}) 
So 
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with the propagator for the neutral mesons 
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Here the a dependence has been reinserted. Note that 
the first term in y ^ p G vv is simply A vv which in the 
case of (3 = 2 exactly cancels the propagator A„„ in (12411 . 
Cancellations like this is a typical property of (3 = 2 
and we shall see several additional examples of terms 
proportional to the factor 8p = 1 — 2/(3, terms which 
only cancel in the case of (3 — 2. Explicitly 
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Taking advantage of the above we find that we can write 
the chiral condensate for general (3 as 

= 1 - I -gN f A vs + 5pA vv 

1 x - {l + S PA ){m 2 +ap 2 )/N c 



V^(p 2 + M 2 v )(p 2 + M 2 vv )TP{p 2 ) 



.(28) 



III. SCHWINGER-DYSON EQUATIONS AND 
VIRASORO CONSTRAINTS 



Having established the equivalence of the replica and 
the supersymmetric methods in all classes of chiral 
symmetry breaking we now turn towards calculating 
Schwinger-Dyson equations that govern the behavior 
of the partition functions in the three classes of chi- 
ral symmetry breaking. That is, we turn from (par- 
tially) quenched chiral perturbation theory to ordinary 
chiral perturbation theory, but the results are also useful 
for the quenched approximation. It turns out that the 
Schwinger-Dyson equations can be expressed in a uni- 
versal way. 

The method is here shortly outlined. For greater de- 
tails see Appendix^] When working with ordinary QCD 
((3 = 2) the finite volumepartition function is an integral 



over the unitary group 



ic partit: 
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d[/(detf/) 1 ' e ^( t ^ t + t/t ^). (29) 
ueu(Nf) 



Here we have included the volume and the chiral conden- 
sate in M. = MVT,o, where M is the usual quark mass 
matrix. 

The partition function is a periodic function in the vac- 
uum angle, 6. This means that the full partition func- 
tion Z(6) may be expressed as a Fourier series, where 
the expansion coefficients, Z„, correspond to the parti- 
tion function in sectors of fixed topological charge, v. 

As discussed by Leutwyler and Smilga |21|. the topo- 
logical charge need not be integer. Letting and tir 
denote the number of left- and right-handed zero-modes 
of the covariant derivative Jp, respectively, the index the- 
orem states that 



1 



n L - n R 



16tt 2 



d 4 xTi(G^G^). (30) 



The trace depends upon the fermionic representation, 
and the trace of the generators, t a , is \i{r)5 ab = Tr[t a t b ]. 
£(r) is the index of the fermionic representation. 
The topological charge is 



1 



32tt 2 



< ^ X {G°n V G a ), 



and thus it is possible to deduce the relation 



l[r) 



(31) 



(32) 



where v is integer. The reason for the appearance of Sp^ 
in eq. 1)32(1 is that for real representations, the reality of 
the Dirac operator leads to all the eigenfunctions of 75 
appearing in pairs. In particular, the number of zero- 
modes is even for (3 = 4. 

For the gauge groups and representations usually con- 
sidered in the literature the relation between v and v 
is simple. For (3=1 one usually considers fermions in 
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the fundamental representation of SU(N C = 2) for which 
£(r) = 1 and v = v. The same is the case for ordinary 
QCD, the fundamental representation of SU(N C = 3). 
However, for (3 = 4 the usual scenario is adjoint (Majo- 
rana) fermions for which l(r) — 2N C and v = N c v (again 
with the gauge group SU(N C )). For simplicity, these are 
also the only representations considered here, but it is 
straightforward to generalize to other representations. 

In the following, the crucial property of the partition 
function, eq. |J2SJ|, is that (det M)~ V Z^= 2 only depends 
on the Nj eigenvalues x iy i = 1, . . . , Nf, of MM^ 0. 
This can be seen from rescaling M. — > V~ X M. and using 
the invariance of the unitary measure under left and right 
multiplication. 

This dependence on the eigenvalues of M. Jvfi only also 
holds in the other classes of chiral symmetry breaking • 
In the case of (3 = 1 the symmetry breaking pattern is 
SU(2Nf) — > Sp(2Nf) and a convenient parametrization 
of the Goldstone fields is UIU* with the partition func- 
tion 

JUeU(2Nf) 

I is the anti-symmetric unit matrix and M. is an arbi- 
trary anti-symmetric complex matrix. Concentrating on 
theories without di-quark terms, we have 



1 

-1 



M 



M 

-M* 



(34) 



where we again define M. = MVSq with M being 
the usual Nf x Nf quark mass matrix. In this case 
(det M)~ u Z@ =1 is a symmetric function of the eigenval- 
ues of MM* 0. 



In the case of f3 = 4 the Goldstone fields can be 
parametrized by UU*, U £ U(N f ). With v = N c v the 
finite volume partition function is 0, 



7/3=4 _ 



dU(det U )X> e i-*«uu*)iM+uu<Mi) ) (35) 
ueu(N f ) 



again with M. = MVEq but this time with M. being a 
symmetric matrix, (det M.)~ P Z^ =4: is a symmetric func- 
tion of M M f 0. 

Useful Schwinger-Dyson equations are obtained in Ap- 
pendix ^ We multiply the integrands in the partition 
functions by Trt a UM^ (J3 = 2), TrfUIlPM^ (J3 = 1) 
and Trt a UU t M^ ((3 = 4), and then subsequently act 
with the differential operator V a . The generators of 
the integration manifold t a are normalized according to 
~ ^u^jk- Since the integral of a total derivative 
vanishes this gives relations among the expectation val- 
ues resulting in the Schwinger-Dyson equations ljA42fl . 
(I A44(l . and (|A46(1 given in the appendix. 



These equations can be written in a remarkably uni- 
form way. Note that the number of generators, Af, of 
the cosets is 1 higher than the dimension of the "origi- 
nal" cosets stated in the introduction because of the pro- 
jection onto sectors of fixed topological charge. For in- 
stance for (3 = 2 the integration goes from an integration 
over SU(Nf) to an integration over U(Nf) by including 
the topological charge. The dimensions of the "original" 
cosets are stated in Table [I] Gathering the differential 
equations deduced in Appendix lAl we see that it is pos- 
sible to write the differential equations universally as 



Nf 



i,3 



IK" 72 K 



n 



-!//2 



1 



(36) 



with v — > v for [3 — 4. This differential equation 
is one of the main results of the present paper, nicely 
combining the seemingly very different integrations over 
cosets into one governing differential equation in terms 
of the fermion masses. All coset dependence is thus in- 
cluded in the single parameter (3. In writing this differ- 
ential equation we have assumed that the mass matrix 
M. is real with positive eigenvalues allowing us to write 

{det M)- = Ui^ v/2 - 



In the following it will be convenient to introduce the 



notation 



M 2 »r r 



(37) 



An alternative way of obtaining useful Schwinger- 
Dyson equations, reminiscent of the one being presented 
here, is to insert a <5^ in the partition function in the form 
of UijUj k , U G U(Nf), subsequently realizing that this 
can be reexpressed as a differential equation, the differ- 
entiation being with respect to the mass matrices. How- 
ever, one should realize that when recognizing (UU'\ as 
a differential operator, in the cases of (3 = 1,4 the fact 
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that the mass matrices are anti-symmetric and symmet- 
ric, respectively, causes complications with regard to the 
diagonal elements. This is especially so for the case of 
(3 = 4, see for instance eq. (IA32|I . This method of ob- 
taining Schwingcr-Dyson equations has for example been 
employed in ref. to obtain the Virasoro constraints 
in the small mass phase. Thus, in eqs. (19) and (28) of 
ref. a factor of 2 is implied for each differentiation 
with respect to a diagonal element. Having deduced the 
differential equation (|36[l . in which all dependence upon 
the mass matrices is through the eigenvalues, we need no 
longer be concerned with this point. 






Coset 


M-l 


1 


SU(2N f )/Sp(2N f ) 


N f (2N f -1)-1 


2 


SU L (N f ) x SU R (N f )/SU(N f ) 


NJ-1 


4 


SU(N f )/SO(N f ) 


N,*F 1 



TABLE I: The number of generators of the cosets in the three 
classes of chiral symmetry breaking. 



A. The small mass phase 

We now proceed to find Virasoro constraints for the fi- 
nite volume partition function by expanding in the eigen- 
values of AiAi^ . For the small mass phase this was first 
done in ref. in the simplest version, namely (3 = 1 
with vanishing topological charge. More recently, Vira- 
soro constraints in the small mass phase have been calcu- 
lated in ref. ■ We present our calculation of the small 
mass Virasoro constraints here mainly for two reasons. 
We wish to demonstrate the ease by which we arrive at 
the Virasoro constraints by using eq. (|36() and to demon- 
strate the correctness of l|36[) by exactly reproducing the 
results of 0. 

In the small mass phase a suitable set of expansion 
parameters for the partition function is defined by 



i 

which transforms the derivatives in H36[) into 

E 

k 

^ \ 1 / -t -. \ k—2 ® 



d 

dxi 



4 fc 1 dt k ' 



(38) 



(39) 



y 1 r k+i-2 9 d 

^ Ak+i * dtkdti' 



k.i 



(40) 



Note that using the invariance of the Haar measure as 
well as permutation symmetry of the Xi, we are allowed 
to remove the sum over i [12j. Inserting the derivatives 
in our master differential equation (|36|l and removing the 
sum over i, we obtain 



V 1 (h - 1 W* 9 + V 1 T k+l 9 9 



k _d_ 

Jk dt k 



| J2k=2 J2l=l V s 0jrx\tk-i-^ 
-i//2 



= 0. 



(41) 



Determining the Virasoro constraints is now a matter of 
reading off the coefficients of x\. This procedure is jus- 
tified provided the coefficients are determined uniquely 
and in a consistent way. Uniqueness is satisfied by re- 
quiring the boundary condition Z@ — 1 if all fermion 
masses vanish. Consistency of the equations obtained by 
simply reading off the coeffecients of x\ is a more subtle 
point, all of the expansion coefficients tk being indepen- 
dent only in the Nf — > oo limit. That this procedure is 
consistent has been demonstrated in ref. 11]. We thus 
find the constraints 



n 



-v/2 



r!3 _ 



0. 



(42) 



with the Virasoro operators L@ 



(3, 



N f + -(n5 fj + v) m 



2 dt 7n dtfi—m 

rn—1 



(43) 



Again we have v — > D for j3 = 4. The normalization is 
chosen such that the operators satisfy the Virasoro 
algebra without central charge 



[Li, = (m-n)L 







m, n > 1. 



(44) 



The virtue of the Virasoro operator formalism is that be- 
cause the operators satisfy the Virasoro algebra, in prin- 
ciple all coefficients of a Taylor expansion in terms of the 
parameters (|38|l can be found from L\ and L^. The Vi- 
rasoro operators l|43|l exactly match those found in ref. 
0. For the solution to the constraints we also refer the 
reader to ref. 0. 



B. The saddle-point approach to the large mass 
phase 

Virasoro constraints in the large mass phase for (3 = 2 
with vanishing topological charge were found in ref. [T(i| 
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and solved in ref. |8(. Below we will present the exten- 
sion of these results to include a non-vanishing topolog- 
ical charge. We will also give the first terms in a mass 
expansion of the (3 — 1, 4 partition functions in the case of 
equal quark masses. Analytical results with equal quark 
masses have been found in both of these symmetry break- 
ing classes 22] . The disadvantage of these results is that 
they depend on determinants proportional in size to Nf 
which quickly makes calculations somewhat cumbersome. 
The results presented here, on the other hand, have a 
very simple Nf dependence. 

The difficulty in the large-mass expansion arises from 
the fact that it is not possible to express the partition 
function in terms of a simple Taylor expansion in the in- 
verse masses. However, as noted in the literature [lfllllj . 
for f3 = 2 this becomes possible after extracting a pref- 
actor from the partition function. In this case the pref- 
actor turns out to be given by the saddle-point approx- 
imation, or the classical contribution, of the partition 
function and thus one would expect this to be the case 
in general. The classical contribution is defined as the 
Nf — > oo limit. We have been able to determine the 
classical contribution to the next to leading order only in 
the case of (3 = 2. For (3 — 2 this is the order needed to 
be able to include the topological charge in the Virasoro 
constraints. In the case of /3 — 1,4 however, we have 
found indications that the next to leading order term in 
the saddle-point approximation is itself given by an in- 
finite mass-expansion. Even so, using the lowest order 
correction in this mass-expansion we have calculated the 
lowest order mass dependence of the partition function 
for [3 = 1 , 4 with equal fermion masses and non- vanishing 
topological charge. Using this result we calculate the one- 
loop correction to the quenched chiral condensate. 

Following |Tot Q] we define the free energy, F, in 
sectors of vanishing topological charge by Z K 



equation 



— c N f F 







Defining in general Z@ — e NfF " we find it useful to in- 
troduce 




£ 



N f J dx. 

dxj 



(3Nf ' Xi — Xj \ dxi 



.(48) 



In refs. 0, 0, 0] the classical solution for (3 — 2 with 
v = was found to be 



Xj)- 



(49) 



1,3 



Ignoring the second derivative and &p , both of which are 
suppressed by an order of Nf , we likewise find the solu- 
tion to (|48|l to lowest order in l/Nf 



£^~ ^£ ln (V^ + V%j)- 



(50) 



From studies of the partition function for small Nf and 
with equal masses, using the analytical expressions of ref. 
[22l | we know that f5H)l is too simple; we expect that Sp 
has to be included in the prefactor. We take care of this 
by expanding the free energy as 



Nf 



(51) 



where is defined by (|50|l . Inserting this expansion 
in (|48|l we now find the following differential equation for 



Gr, 



F„ 



(45) 



-v/2 



since then 
operators this means 

d 



)Z% = e N ' G ". For the differential 



e N ^ = N f e N f G -^G v 

OXi OXi 



d 2 



dx 2 



N f e 



N } G V 



N f 



fdG v \ 2 d 2 G v 



V dxi 



dx 2 



(46) 



.(47) 



In what follows we will use the somewhat sloppy termi- 
nology that G v is the free energy. 

The free energy, being an evaluation around the Nf — > 
00 saddle-point, can be written as a power expansion 
in l/Nf. When doing the counting of orders of Nf the 
action counts as order Nf. Explicitly including this fac- 
tor to ease the counting of powers, when inserting the 
free energy into the differential equation (|36|l we find the 



Sp + 2v 



y/XiXj 



4 (v^+v^J) 2 



(3N S 4^ 



dxi (3Nf t— 1 yfxl + ^/xj dxi 



£ 



XiXj ( d d 
[3Nf ^— ' Xi — Xj \ dxi dxj 

-.(1) 



G {1 



(52) 



where we note that G v vanishes for (3 = 2 with vanishing 
topological charge since 62 = 0. 

At this point a few comments are in order. All terms 
in (|52|l have to cancel order by order in x, x £ Xi, i = 

1, . . . , Nf. Thus diG^ only has terms of order less than 
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or equal to 0(x 1 ). Otherwise cancellations would only 

happen if did 1 ^ had terms of all positive orders of x. 
From this consideration it is clear that the term of or- 
der y/x in ijol'l) should be cancelled by the term of or- 
der x 3 / 2 diG^\ This provides a condition for Gv which 
should then be inserted in l|52[) . If this expression fails 

to cancel the remaining terms we conclude that diG^P 
also has an x~ 3 / 2 part and the process starts over again. 
Hopefully the process terminates after a finite number of 
iterations. 

From these considerations we conclude that 



d 5p + 2v 



— G« = - 
dxi v 

G (1) = - 



Ax, 
Sp + 2v 



(53) 



where in the last step we used the fact that G^ should 
be a symmetric function of the Xi. Unfortunately the 
iteration process does not stop within the first three lev- 
els of iteration. Indeed, the results become increasingly 
complicated - thus we choose to stop the iteration at this 
point, using H53[l as our only correction to l|50(l . It turns 
out that this is enough to calculate the partition func- 
tion if one requires the fermions to have equal masses. 
However, since all remaining terms are proportional to 
5p, is precisely sufficient to calculate Virasoro con- 
straints in the (3 = 2 case with non- vanishing topological 
charge. Furthermore, the remaining terms are indepen- 
dent of v. Thus, in all classes of symmetry breaking, to 
this order the v dependence of the free energy is simply 
given by eq. 

First we consider (3 — 2 with non- vanishing topological 
charge. We write the partition function as 



h3 
-u/2 



J\x- v/2 Y{{t k }). (54) 



In the second line we have removed the explicit in- 
dependence of the action, and returned to our original 
normalization. We define the set of variables relevant to 
a large mass expansion differently from the small-mass 
expansion 



tk 



2k 

22fe+l 



-Tr ({MM 



U-k-1/2 



2k 



1? 



-fc-1/2 



(55) 



with the normalization chosen for convenience. Inserting 
(|54|l in the governing differential equation, eq. i|36|) an d 
reading off the coefficients of powers of Xi as in the small 



mass expansion we find the Virasoro constraints {y — > v 
for (3 = 4) 

L P n = 2 Y{{t k }) = 0, n>0, (56) 
with the Virasoro operators 



£fc 2 = 5> + l/2)f ( 

fc=0 

d 

din 



d 



dth+n 

n > 1. 



1G 



if 



4 ^ dt k -idtn-k 



(58) 



Again these operators satisfy the Virasoro algebra with- 
out central charge. Uniqueness of the solution is ensured 
by the boundary condition Y({t k }) = 1 for all expan- 
sion coefficients vanishing (the limit of infinite fermion 
masses). Comparing these Virasoro operators with the 
Virasoro operators obtained in ref. 01 we see that 
the inclusion of a non- vanishing topological charge only 

Q— 2 

changes the Virasoro operator Lq . However, this still 
has profound consequences for the power expansion of 
the partition function, as can be readily checked. 

To determine the first corrections due to including a 
non-vanishing topological charge, we proceed by expand- 
ing the function Y as 

oo 

Y({t k }) = 1 + J2 c ktk+ e *i,fa t *i*fca + -"-(59) 

fc=0 0<fci<fe 2 

The coefficients are easily obtained from the Virasoro 
constraints and we find to fourth order in the inverse 
masses 



i,j v 

l + l^i^ Tr{MMt) -l/2 



(l-4i/ 2 )(9-4^ 2 ) 



128 



(Tr(MM t )~ 1/2 ) 2 

2VT7 _ A,,2\ 



(l-4, 2 )(9-4, 2 )(17-4, 2 ) 

3-1024 V V ' ' 



(l-4i/ 2 )(9-4^ 2 ) 



Tr(MM r ) 



f\-3/2 



3 • 128 

(1 - 4i/ 2 )(9 - 4t/ 2 )(17 - 4;/ 2 )(25 - Av 2 ) 

3 • 32768 
x(Tr(WM f )- 1/2 ) 4 
(l-4^ 2 )(9~4i/ 2 )(25-4^ 2 ) 
3 • 1024 

xTr(XX t )- 1/2 Tr(A!Al t )" 3/2 + 



(60) 



The partition function for (3 = 2 is known in closed form 
to be IH 



z?r 2 



detA{{^}) 
A({x}) ■ 



(61) 
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A is an Nf x Nt matrix defined by 

- xt 1)/2 Iu +j -i(V^), (62) 

with /„ a modified Bessel function. A is the Vander- 
monde determinant 

N f 

A({x}) = Ylfa-Xj) 

i>j 

= detKxiY- 1 ). 



(63) 



Expanding l|61l) for Nf = 2, which is sufficient to reveal 
most subtleties, we find complete agreement between this 
closed expression and our expansion, equation (|60fl . 

Turning now to the other classes of chiral symmetry 
breaking, we find the partition function for equal fermion 
masses. This is performed as before, by writing Z@ as 
the prefactor times a power expansion in the masses, that 
is 

1,3 1 ■> k 

Inserting this expansion in (|36|l and setting all fermion 
masses equal, after some algebra, we obtain 

-2c, = l + (*/-l)2W-4^ (65) 
8 

where, as expected, we find that the partition functions 
are equal for Nf = 1, a feature which should be present 
not only to this order in the mass expansion but to all 
orders. We have here returned to our original normal- 
ization, removing the explicit Nf dependence from the 
partition function. Seeing how the v dependence of the 
free energy was found to be very simple to the order 
examined, it should not come as a complete surprise 
that the ^-dependence is independent of 8p, apart from 
the difference in prefactors. Explicitly writing the parti- 
tion function for equal masses to lowest order (defining 

X = Xx = • • • = Xjsr f ) 

Z I3 = e N f y/S x -N}/20 x -N f Sp/4 

^ i+ l + ( iV / -l W -4,^ + ^V 6) 



where we as usual let v — > v for f3 = 4. 

This result is just what we need to calculate the fi- 
nite volume quenched chiral condensate in the large mass 
limit. To this end we replace Nf with N v or, formally, 
extend the partition function with one replica set of N v 
fermions and then let Nf — > 0. Changing variables 
to fx v = y/x (with the equal fermion masses given by 
fEoV) eq. can also be written 



%1 = hm i-^m^. 

E w„-o N v d^ v 



We thus calculate the quenched chiral condensate in sec- 
tors of topological charge v (v — ► v for (3 = 4) and in all 
classes of chiral symmetry breaking 



So 



= 1 - 

quenched 



2fi v 



1 - 2^//? - 4j/ 2 
8/4 



(68) 



While the quenched chiral condensate has been previ- 
ously calculated in the case of f3 = 2 [17|. interestingly 
for (3 = 1,4 our results show that we have to include a 
term of order ©(/it" 1 ) in these classes. This feature turns 
out to be very important in the matching of the chiral 
condensate in the two finite volume perturbation schemes 
considered in the next section. 



IV. QUENCHED FINITE VOLUME CHIRAL 
CONDENSATES 



So far everything we have calculated has been in the 
large, albeit finite, volume, in which the usual momen- 
tum expansion is applicable. Now we turn to the other 
finite volume regime, in which the correlation lengths of 
the Goldstone modes are much larger than the volume of 
the box. In this regime the usual p-expansion of chiral 
perturbation theory breaks down due to the propagation 
of zero-momentum Goldstone bosons. Thus another ex- 
pansion scheme, known as the e-expansion |l5i l24l |. is 
required. In ref. [17| it was shown, in the case of j3 = 2, 
that a region exists in which the two expansion schemes 
overlap. The main result of this section is that matching 
is not fortuitous, the result carries over to the (3 = 1,4 
cases. We demonstrate this by calculating the quenched 
chiral condensate in the e-expansion and comparing with 
the corresponding result l|28|) of the p-expansion. Thus 
we conclude that matching of the two expansion schemes 
is possible in all classes of chiral symmetry breaking. 

Again we use the 0{p 2 ) Lagrangian which we param- 
eterize as in (J3J), i.e. 



£ eff ~ 



^Tr(c) M C/ta M C/) - ^-Tr(M^(U + tft)) 



ra, 
~2N, 



2N, 



(69) 



(67) 



The virtue of the e-expansion is that it counts zero- 
momentum modes as order 0(1) while non-zero modes 
count as higher orders in e. Letting the four- volume be 
V = L 3 /T and treating the fermion masses as small in 
comparison with T and 1/L one fixes the counting of or- 
ders of e by setting 1/L = 0(e). The fermion masses rrii 
are then counted as 0(e 4 ) which implies Mjj ~ 0(e 2 ). 
Thus all graphs which exclusively involve zero momen- 
tum propagators are counted as being of order 0(e°) 
|L5| . see for instance eq. i|10|) . In the unquenched the- 
ory, adding a mass term for the non-zero modes provides 
a Gaussian damping factor for these modes which en- 
sures that the non-zero modes count as 0{e 1 ) |15| . In 
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the quenched theory the counting is not as simple since 
the mass term m Q is a free parameter. But following 
[T?l l2fij we perform a simultaneous expansion in e and 
TOq/ (N c F 2 ), as a way of avoiding having to evaluate all 
orders of mo for each order of e. 

The idea, adopted from ref. j is to collect the zero- 
momentum modes in the constant matrix Uq € G/i? and 
write 



C7(x) 



(70) 



where the field £(x) G G/if contains the non-zero mo- 
mentum modes. For the non-zero modes we choose the 
same parameterization as in (JSJ, for j3 = 1, and l|17|) . for 
/3 = 4. That is (/3 = 1) 



1 



71 U f T 



(71) 



with (/) Hermitian and ^ anti-symmetric complex fields, 
and (J3 = 4) 



An j=A 12 



^A hNf+Nv \ 



A 



,(72) 



Nf+N v ,N f +N v J 



with £ being a real, symmetric matrix. 

We now wish to perform this change of variables and 
integrate over the zero modes and the non-zero momen- 
tum modes separately. Thus we have to perform the 
change of variables l|70[l in the partition function. The 
Jacobian of this change of variables is, to lowest order, 
just a constant, corresponding to a shift of the vacuum 
energy. Noting that integrals of odd powers of £ van- 
ish, with this change of variables the lowest order action 
becomes 



S = / <T xTr 



-d^{x)d^{x) 



+ ^o + ^)f + ^E(x)d^(x) 



~2~ 



M{U G + Ul) (v-js Jdxi 2 



, (73) 



where Uq = e l ^ 2ipo l Fl3 , and So = Trt^o- The mass matri- 
ces M are defined as in eqs. © and i|T5)l. 

Since we are once again interested in calculating the 
one loop correction to the chiral condensate we integrate 
out the non-zero momentum modes, £(x). To this end 
we define the non-zero momentum modes 



A, 



V ^ 



1 



^ = ^E G ( M S)- 



(74) 
(75) 



The two-point functions for the off-diagonal mesons are 
related to the full propagators by 



1 



(76) 



The propagator of the diagonal mesons, G(M^), is de- 
fined as in eq. (|25|l . 

Thus integrating out the fluctuating fields we find 
the fully quenched one loop correction to the effective 
Lagrangian for the zero momentum modes by putting 

jvy-o, 



lim — 

w„->o N v 



2Fj 



Tr 



M{U + UI)[ dx£(x) 



v 1 VE 
lim — 

n v ^o N v 2F 2 



Tr(M(U + U^)) 



= lim 



1 VE 



1)A, : 



G v 



n v ^o N v 2F 2 



Tr(M(Jf + Ufo) 



x ( S () A VV - (1 + 5[3 A )—A VV + 

•(1 + 6p,i)(rn% - aM* v )d M 2 A vv ) .(77) 



In the last equality we have made use of the fact that the 
first part of the diagonal propagator is the same as for 
the off-diagonal propagator. That is 



Gij — Aij 



(6 0}1 + l)(ml + ap 2 )/N c 



^ (.P 2 +M 2 l )( P 2 +M 2 ] )TP{p 2 y 



(78) 



Equation (|77|l is very useful, since from the correction to 
the zero'th order approximation we see that the one-loop 
improved effective partition function can be written in 
the same way as the zero'th order approximation 



70 _ 



G/H 



dUoe hv^TrM(u +u^-^fs^ (79) 



by simply changing the mass eigenvalues into 
1 



{N C 5 A VV - (1 + 5^)aA v 



F 2 N C 

+ (1 + fy,i)(m§ - aM* v )d M 2 v A vv } 



(80) 



Notice that we again define fi v = m v VT,Q. 

Thus we can exploit the usual formula, Q , to find the 
one loop improved quenched chiral condensate 



£(Mu) 

So 



1 dfi' v d 



lim 



(81) 



Now we make use of our calculation of the chiral con- 
densate in sectors of fixed topological charge. When cal- 
culating the chiral condensate, we have to sum over topol- 
ogy |26j . To perform this summation we follow the idea 
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of ref . [lTj and include an arbitrary vacuum angle in the 
partition function by absorbing it into the field S . Next 
we perform a Fourier transform projecting the vacuum 
angle onto sectors of topological charge. 

Identification of the topological susceptibility |2(j is 
performed in Appendix iBl The result is 



/ 2\ /, r ,m 2 VF 2 



2N r . 



The summation over topologies is 



oo 

£ 

u— — oo 



Z^)Eg(^) 
ZOO*) ' 



(82) 



(83) 



Inserting the quenched chiral condensate in sectors of 
topological charge calculated in eq. JBSJl, to the order 
considered we find 



= 1 - 



So 

h_ 

2fj, v 



l-26p/(3-A( V 2 ) 
8 M 2 



(84) 



Now we look for a size of the volume V in which the 
two finite volume regimes match [13. Obviously JH1 
cannot be a simple perturbation of the infinite volume 
theory, since \i v = m v VT,Q in this regime is finite while 
[v 2 ) ~ V . Thus the requirement that the last term in 
is a small perturbation is 



v » (i + ^,0 



F 2 ml 
AN c mlY,l 



(85) 



The requirement of this extreme finite volume regime is 
that V < M~ v 4 = F 4 /4m^I]2. Both limits can be met by 
letting N c — > oo while keeping all other quantities fixed, 
since then F 2 > (1 + 5 f 3.i)ml/N c . 

Having established a region in which a match is possi- 
ble, let us first review the result i|28|) for the chiral con- 
densate in the large, but finite, volume. Rewriting this 
for the quenched limit we obtain 



So 



1 



= I - — — {NcSpAw - (1 + 8f3.i)aA vv 

aM 2 v )d M;v A vv } .(86) 



F 2 N C 
+(l + ^,i)(m§ 

In the extreme finite volume we also have all the in- 
gredients necessary for calculating the quenched chiral 
condensate. Inserting eqs. (|8(J|I and (|84|l in eq. (|81|l we 
find to leading order 



- 1- 



00 



VF 2 Ml 
1 

" F 2 N r 



+ (i + <W) 

{N c 5pA v 



N C F 2 M$ V V 



1 - 



-(1 + fy,i)(mg - aM 2 v ) 9m$ v A vv } 
■ {N c 5pA vv - (1 + ^,i)«A OT 



(1 + 8 PA )(m 2 - aM 2 v )8 M 2 v A vv } .(87) 



In this calculation, we have made use of the fact that 
(v 2 ) > 1 while [i v - O(l). 

Comparing the results l|87|) for the extreme finite vol- 
ume and (|86() for the ordinary finite volume, we find the 
promised overlap of the two perturbative schemes. This 
is an important consistency check on our calculations and 
it is a very attractive feature of this calculation that it is 
precisely the term proportional to l/fj, v in the quenched 
chiral condensate, absent in the (3 — 2 case, which makes 
this matching possible. 



V. CONCLUSIONS 

We have analyzed the three patterns of chiral symme- 
try breaking relevant to chiral perturbation theory, in 
particular (3 = 1 and (3 = 4 corresponding to fermions in 
a pseudo-real and real representation, respectively. 

In the usual p-expansion of chiral perturbation theory, 
we have demonstrated how the replica method is equiv- 
alent to the supersymmctric method by comparing the 
Feynman rules to one loop order for the (3 = 1,4 classes 
of chiral symmetry breaking, and seeing how the calcula- 
tions proceed it is fairly obvious how this generalizes to 
higher order loop calculations. 

Next, using Schwinger-Dyson equations, we have ob- 
tained a governing Schwinger-Dyson equation for the par- 
tition function in all classes of chiral symmetry break- 
ing and in sectors of fixed topological charge. We have 
demonstrated how it is possible to write this differential 
equation in terms of the mass eigenvalues only, with all 
dependence upon the class of chiral symmetry breaking 
being through the Dyson index, (3. 

We have utilized the governing differential equation to 
determine Virasoro constraints in the small mass expan- 
sion, which exactly matches those derived in ref. 0- 
We have also calculated Virasoro constraints in the large 
mass expansion of the (3 = 2 partition function in sec- 
tors of fixed topological charge. We note how the inclu- 
sion of the topological charge in the Virasoro constraints 
becomes possible by simply expanding the saddle-point 
contribution of the free energy to next to leading order 
in l/Nf instead of only keeping leading order terms. 

We find indications of how the next to leading order of 
the saddle-point approximation for (3 = 1, 4 seems to con- 
sist of an infinite series in the fermion masses. However, 
we note that it is still possible to find the first terms in a 
power expansion of the partition function in these classes 
of chiral symmetry breaking if one puts the masses equal. 
Having thus obtained the equal-mass partition function 
in sectors of fixed topological charge to the order needed, 
we calculate the quenched chiral condensate. We find 
that the mass-dependence of the quenched chiral conden- 
sate is different for real and pseudo-real representations 
of the fermions compared with the formerly known result 
for fermions in complex representations. 

In the other relevant finite volume, the volume smaller 
than the correlation lengths of the Goldstone bosons, we 
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utilize the e-expansion to calculate the first correction to 
the quenched chiral condensate. It turns out that it is ex- 
actly the above mentioned difference in mass-dependence 
of (3 = 1,4 which makes the quenched chiral condensate 
of this finite volume match the corresponding quenched 
chiral condensate of the ordinary p-expansion. Thus we 
find a region in which the two perturbative expansions 
of chiral perturbation theory, the e-expansion and the 
p-expansion, precisely match. 
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APPENDIX A: SCHWINGER-DYSON 
EQUATIONS 

In this appendix we will find Schwinger-Dyson equa- 
tions in all three classes of chiral symmetry breaking. 
Theories with complex fermions. such as ordinary QCD, 
have been considered before 14], but since the calcula- 
tions here are basically a less complicated version of the 
calculations in the (3 = 1,4 cases, we will line them up. 
Also, this will be done in a somewhat different way from 
the literature. 

To avoid an unnecessarily complicated notation we 
consider first theories with vanishing topological charge. 
The topological charge is easily included at the end of 
the calculations. 



these is that X/ a UW — 0. An explicit representation of 
V a is 



(A3) 



We normalize the generators of the unitary group accord- 
ing to 

1 



ij L kl 



il°jk- 



(A4) 



Getting a useful Schwinger-Dyson equation is now a 
matter of choosing the function F(U) in a clever way. It 
turns out that with the choice of F(U) — UAi^ we arrive 
at the equation considered in 01 

Nf (TrM^U) = i (TiMM 1 - TrUM^UM^) , (A5) 

where we have used i|A4|l . 

Now we follow the idea of ref . [14| . We know that 



(TrM^U) = 2 / M\ 



d 



l3 dMl 



2 <£>!>, (A6) 



(TrUM*UM*) = 4(Mt.A4t— ^-r- 



4(D 2 



(A7) 



To derive how the differential operators D\ and D 2 act 
on the partition function in terms of the eigenvalues we 
use the fact that Z^ =2 only depends on the Nf functions 



cj) p = Tr((MM^) p ) 



(A8) 



1. (3 = 2 

We start out with the partition function in sectors of 
vanishing topological charge 



= 2 



/ dUe^ T < MU " +M " u ). (Al) 



As discussed in section Mil the crucial property of this 
partition function is that it only depends on the Nf eigen- 
values Xi of AiAi^ . To arrive at the Schwinger-Dyson 
equations, we utilize the fact that the integral of a total 
derivative vanishes. Thus 

= / dUTrt a V a (F(U)ei( Mui+Miu )) 

Trt a V a F(U) 
1 



Tr(i a F(J7))V a Tr (MU^ + M ] U)\ . (A2) 

The left derivative V a on U(Nf) is defined through the 
relations V a U = it a U and V a U^ = -illH a where t a are 
the generators of U(Nf), A simple consistency check on 



p = l,...,Nf. This follows from the fact that Z@ =2 de- 
pends only on the eigenvalues as well as on permutation 
symmetry among these. Using the chain rule we find 



D ^ =2 = E^^ 



(A9) 



so we only need to know how D\ operates on the basis 
functions, i.e. 

D^ P = pMljdMM^M)^ (A10) 
= Wv (All) 

Thus we see that 

d 



Di = x t 



dx,. 



In the same way 



D 2 <f> p = ^2 p(f> p >(f>p-p> 
P '=i 



E 



(A12) 



^Mpxf-'+K' 1 )' ( A13 ) 
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and following [lij we see that the second derivative is 
equivalent to 



group we can follow the same approach as in the case of 
(3 = 2. We know that 



Do 



E 



a 



_d_ 



(A14) 



As noted in 0] this procedure needs further justification. 
From the invariance properties of the partition function 
we know that Z^ =2 is a function of the basis functions 
only. Accordingly we should be careful when concluding 
which form the double differential operator takes in terms 
of the eigenvalues. Expanding the partition function in 
terms of the basis functions we see that in principle things 
could become complicated from the chain rule when D 2 
acts on a product of the basis functions. We can check 
that this is not the case by letting our "guess" (|A14|) act 
on a product 4> p 4> q and comparing this with the result 
obtained by letting the known correct matrix form of the 
derivative l|A7jl act on 4> p 4> q . Since we are dealing with a 
double derivative no new complications arise if we look 
at a term (j} p (j) q <j) r ... in the expansion of the partition 
function. This test, in fact, gives the correct result. 

Using (|A5|) we collect the terms and find the following 
Schwinger-Dyson equation for QCD with fermions in a 
complex representation: 



\dl + NfXidi + 



X 1 X- j 
Xi X- 



(ft - d 3 ) 



iE^ 2 - 



rP=2 



(A15) 



0= / dUTrt a V a (F{U)e^< MW+jCiv ')) , (A17) 

and if we, inspired by the (3 = 2 calculation, choose 
F(U) = VM) this results in 

(2N f - l)TrMtyJ> = l-(Ti-MM^ 

-TcVM*VM*y (A18) 

It is important to realize that since both V and JVv are 
antisymmetric, this equation should be treated a bit dif- 
ferently from the (3 = 2 case. Being careful, we define 
the differential operators 




d 



(A19) 



(TrVM^VM^ 



4/ J2 MjMl 



i,j,k,l 



d d 



4(D 2 



(A20) 



Again we let these act on the basis functions <j> p ,p = 
1, . . . , Nf. The basis functions are defined in the same 
way as in the (3 = 2 case but differ, since the matrices in 
this case are 2Nf x 2Nf. Asuming that the theory has 
no di-quark terms, the mass matrix has the form 



M = 



M 
-M* 



(A21) 



In the case of fermions in a pseudo-real repre- 
sentation the Goldstone bosons live in the coset 
SU(2Nf) I Sp(2Nf). This coset can be parametrized by 
{/J/7', where / is the antisymmetric 2Nf x 2Nf unit ma- 
trix 



1 

-1 



and U € U(2Nf). The partition function with vanishing 
topological charge is 



Z' 3 



=1 



dUe kTr{M^UIU<+M(UIUy)^ (Alg) 
UeU(2N f ) 



where the mass matrix M is an arbitrary complex an- 
tisymmetric matrix. Note also that the combination 
UIU* = V is antisymmetric, which will turn out to be 
very useful in the following. Again, the partition func- 
tion is a function of the eigenvalues of .MA^ and since 
the Goldstone modes are parameterized by the unitary 



which results in the basis functions (again letting Xi, i 
1, . . . , Nf be the eigenvalues of the matrix M.M.') 



4> P = Tv{{MM^) p ) 

= 2E*?- 



(A22) 



Though the sums in eqs. I|A19|) and (|A20|) do not run over 
all indices, when doing the actual derivatives the results 
turn out rather nicely. Indeed for the first derivative we 
find 



d 



t 

ij dMl 



j2m 



-((MMl) p - 1 M) i 



2p<() p 



(A23) 
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from which we recognize 



£>! = 2 Xl — . 

oxi 

Turning to the second order differential operator 



(A24) 



i,j,k,l 



dMl 3 dM a 



P '=i 



Mj 



-(x: 



p-i _ -p- 1 ^ 
x i / 



+8p(p-l)X; a! ?-4p(p-l)X; !E ?, (A25) 



from which we conclude that 

3j ^ <JT ^ 



£> 2 



2a; 2 92 , iy ggj ( d (A26) 



Mi 



Again we have to verify this double derivative on a prod- 
uct ifipfiq- The test produces exactly the necessary result. 

Collecting the terms using eq. I|A18|I we find the dif- 
ferential equation 



Mi J 



3. (3 = 4 



(A27) 



For fermions in a real representation the coset is 
SU(Nf)/ SO(Nf). Once again this can be parameter- 
ized by elements from the unitary group, in this case by 
UU . Here the mass matrix is an arbitrary symmetric 
complex matrix and the par tition function with vanish- 
ing topological charge is [2l], |22| 



is that in this case the diagonal elements of the mass 
matrix are not equal to zero. Again we carefully define 
the differential operators 



Mi 



= (Di) , 

\ fjt,l 9M kj dM il 

MW= 
frf dM], dMl 



(A31) 



= (D 2 



(A32) 



The partition function only depends on the eigenvalues 
of A4A4^ and since we have returned to working with 
Nj x Nf matrices the <p p , p — 1, . . . , Nf are exactly the 
same as for = 2. That is 



In this case we find that 

D\(j)p = 2pcj) p , 



and thus 



D 1 = 2xi 



d 



(A33) 



(A34) 



(A35) 



The second derivative looks complicated but after some 
algebra we end up with the nice expression 



z p=i = / dUe^( Miuut+M ^ uu ^).(A28) 

JU£U{Nf) 

Proceeding in the familiar manner, with W being the 
symmetric matrix UU , we find 

= / dUTrt a X7 a ( F (U)e^ M ' w+Mwr )\ . (A29) 

Defining F(U) = WAi^ we find the Schwinger-Dyson 
equation 

((N f + ljTiM^W) = (TiMM^ - TiWM^WM^)(A30) 

Having seen the techniques in the anti-symmetric case 
of = 1 this is much the same. The only complication 



D 2 <t> p = 2p ^ 4> P ,( t>p-p' + 2p(p~ 1)0 P , (A36) 
p'=i 

from which we deduce that 

^^Strt)^ (A37 » 



dxf 



Mi 



Now we are able to collect the terms and we end up with 
the result 



Mi 



(A38) 
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4. Including the topological charge 



the differential equation 



Including the topological charge in the Schwinger- 
Dyson equations is now quite simple. In all classes of 
chiral symmetry breaking the topological charge enters 
the partition function as a change of measure. Again we 
make use of the fact that the integral of a total derivative 
vanishes. Thus we get an extra term in the Schwinger- 
Dyson equations corresponding to the derivative of the 
change in measure. 

More specifically, consider first (3 — 2. Allowing a non- 
vanishing topological charge changes the measure in (|A1|) 
by dU — > dU (det U) v . Thus we need to know the deriva- 
tive of det U 



V a (detC/) = i{t a U) 



d 



'',r 



detU 



dUij 

iieU^jUj^ det U 
iTr(i Q ) det U 



(A39) 



since 



det(U + SU) = det(f/(l + U^SU)) 

= dct[/(l + Tr([/- 1 (5Lf) + ...).(A40) 



Thus we see that the expectation value equation (|A5|) 
changes into 

(N f + u)( TrM ^U) = ^( TrM M f -Tr U M f U M f ) , ( A4 1 ) 

which only alters the coefficient in front of the differential 
operator D±. Modifying the partition function, of course 
also implies that the dependencies are changed. As dis- 
cussed in section ITTT1 (det M)~" Z^~ 2 is a function of the 
eigenvalues of MM^ only, and all results of the last sec- 
tion are thus valid in sectors of non- vanishing topological 
charge provided we change the coefficient in front of D\ 
and let the differential operators act on (det M)~ v Z@ =2 
instead of just Z^ =2 . Thus we find 



\d\ + (N f + v) Xi d i + £ (a, - a 



in 



(det M)- U Z^= 2 = i^x !; (detX)^Z^= 2 .(A42) 



Turning to (3 = 1 where the measure also changes by 
dU — ► dU {det U) v we find that the expectation value 
equation i|A18(l changes into 



(2JV> - 1 + v)TrM.W 



\ ( TtMM^ - TrVM^VM^ 



(A43) 



In this case we need to replace Z^~ x with 
{detM)- v/2 ZP =1 = (det M)-" 'Z^ =1 and we find 



■<■:»: + pN f - 1 + uirji + 2E„ 3 i^r- (A - dj) 

x (det M)- V ZP =1 = ^i^i^MyZ^ 1 . (A44) 



Finally, for (3 = 4 the change of the measure is dU - 
dU(detU) 2p , with v = N c v. Hence eq. (|A30|) becomes 



= %(TrMMl -TrWMtWMl). (A45) 
Z /3=4: should be replaced by (det M)~ U Z^~ 4 and thus 



x{detM)- p zP- 4 = i^ ia;i (detM)- p Z^- 4 .(A46) 



APPENDIX B: THE TOPOLOGICAL 
SUSCEPTIBILITY 

In this appendix we calculate the topological suscepti- 
bility. This we calculate by allowing a non- vanishing vac- 
uum angle 6 and subsequently projecting onto sectors of 
fixed topological charge by means of a Fourier transform. 
We consider the classes of symmetry breaking separately. 

The case of f3 = 2 has been considered previously [13 
by means of the same method as we will use here. The 
result found was 



1 0=2 



mlVF 2 
2N C 



(Bl) 



Matching of the two finite volume regimes requires that 
the dependence of the topological susceptibility upon the 
Dyson index follows the same pattern as the last (non- 
trivial) part of the diagonal propagator. Thus we expect 
the topological susceptibility to be twice as large as for- 
mula 1|B1(1 in the (3=1 case while being identically the 
same in the (3 = 4 case. 



1. (3 = 1 

A non-vanishing vacuum angle can be absorbed into 
the mass-matrix M by letting M -> e i6/N fM plll^. 
We can include the vacuum angle in a shifted So- To see 
this consider how the change of the mass matrix affects 
the combination 



ReTrM — > ReTr 



Me ie/N} e -' l ^/ Fl g-'v^o/Fi 



(B2) 



from which we find that we can include the vacuum angle 
in a change of variables 



0Fi 
V2N f ' 



(B3) 
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Taking the trace of this equation we obtain (note that Finally we perform the Fourier transform (for large 
the dimension of the coset is 2Nf x 2Nf) volumes) 



S' = H - V29Fl (B4) 



Z » ~ 2n 



dde- lvB I dU exp 

lsU{2N f )/Sp(2N f ) 



Nr. 



2m$VF{ 



-v 2 N /4nigVF? 



U(2N f )/Sp(2N f ) 



dU (detU )- v/2 



exp 



-TrM(U + W 



(B5) 



from which we identify the topological susceptibility 
2mlVF? _ mlVF 2 



<- 2 > 



Nr. 



N c 



(B6) 



As expected, this is exactly twice the result from the 
(3 = 2 case. 



2. (3 = 4 

The (3 — 4 calculation proceeds in the same way. In 
this case the vacuum angle can be included in the mass 
matrix by letting M -» e l0/N " N fM [Ul^. A S ain this 



can be included in the field of constant modes, this time 
by the substitution 



<^o -» Vo 



6F 4 



V2N c N f ' 

Thus for the traced equation we find 

6F 4 
V2N C ' 



Hq — > S = H 



(B7) 



(B8) 



The projection onto sectors of fixed topological charge is 
performed as before 



7/3 = 4: 



— I dde- w6 
2tt 



dUo exp 



SU(N f )/SO(N f ) 



1~ U m V ^2 



-TtM(U + U q 



2N r 



oc 



-* 2 /2(* 2 ) 



df/ (detJ7 )" P exp 



U(N f )/SO(N f ) 



-TrM(U + Ul) 



(B9) 



with the topological susceptibility given by Again we define v = N c v. This result is also as expected. 



^ I ~ ^2NT ~ ^2NT- (B10) 
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